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Abstract 

We study the expectation of linear eigenvalue statistics of matrix models with any 
(3 > 0, assuming that the potential V is a real analytic function and that the corresponding 
equilibrium measure has a one- interval support. We obtain the first order (with respect to 
n^ 1 ) correction terms for the expectation and apply this result to prove bulk universality 
for real symmetric and symplectic matrix models with the same V . 

1 Introduction and main results 

We consider ensembles of n x n real symmetric, hermitian or symplectic matrices M with the 
probability distribution 

P n (M)dM = Z~l exp{-^Try (M)}dM, (1.1) 

where /3 = 1,2,4 corresponds to real symmetric, Hermitian, and symplectic case respectively, 
Z n a is a normalization constant, V : M — > M + is a Holder function satisfying the condition 

V(X) > 2(1 + e) log(l + |A|). (1.2) 

The joint eigenvalue distribution which corresponds to (II. ID has the form (see |12j ) 

n 

PnA^-^n) = Q-^U e ~ nmXl)/2 II l A *" A /> (!'3) 
i=l l<i<j<n 



where 

n 

-n/3V(\ l )/2 
i=l X<i<j<n 
This distribution can be considered for any /3 > 0. We denote 



,. n 

Qn, P = j n 



- [] \Xi-\jfdXL..d\n. (1.4) 



%{(• • • )} = j (■ ■ ■ K/3(Ai, A n )dAi, . . . d\ n , (1.5) 

and 

Pi$(\i, Xi) = / p„ i/ 8(Ai,...Aj[,Aj[ + i,...,A n )rfAj[ + i...dA n . (1.6) 
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It is known (see [21 [10]) that if V' is a Holder function, then the first marginal density p[ n ^{X) 
converges weakly to the density p(X) (equilibrium density) with a compact support a. The 
support a and the density p are uniquely defined by the conditions: 



v(X) := 2 J log |ju — X\p(p)dpL- V(X) = supv(X), X G a ^ ^ 

v(X) < supu(A), A a, a = supp{p}. 

If we consider the linear eigenvalue statistics of a smooth test function / 

n 

M,[/] = X)/(A0, (1-8) 

i=i 

then the above results of [21 [TO] mean that 

lim E /3 {n-W n [/]}= lim / f(X)p\ n h\)dX = [ f(X)p(X)dX, 
lim {|n-Wn[/] - E^n" 1 A^[/]}| 2 } = 0. 

n— >oo 

Moreover, in [2] some rather rough bounds on the rate of convergence were found 



f(X)(p^(X)- P (X))dX 



<C\\f\\l /2 ||/'||fn-V2logVV (1.9) 



E^jln-^m-E^n- 1 ^!/]}! 2 } < 1 1 / 1 1 2 1 1 / r 1 1 srz" 1 logn. 

Here and below we denote by 1 1 . 1 1 2 a standard I? (c e )- norm, with cj £ being the e-neighborhood 
of the support a with sufficiently small e. 

In the case of j3 = 2 these bounds can be improved considerably. It is a simple exercise 
(see e.g. [TO]) to show that for any V satisfying (j 1 . 2 1) (not necessary Lipshitz) the l.h.s. 
of the second inequality is 0(n~ 2 ), but for other /3 this fact is not proven yet. With the 
first inequality of (|1.9p the situation is similar. It follows from the results of [3] that for 
real analytic V the l.h.s. of the first inequality of fjl .9|) is 0(n~ l ) (see also pQ where the 
asymptotic expansion with respect to was constructed in the case of even real analytic 
V and one or two interval support a). Unfortunately, similar results are not found for j3 7^ 2 
in the general case of a till now. 

Bounds of the type fjl .9|) are interesting not only themselves. They have a lot of very 
important applications, which includes Central Limit Theorem (CLT) for linear eigenvalue 
statistics, the asymptotic for logQn^, etc. One of the most important and interesting appli- 
cations is that to the universality problem for (3 = 1,4. Universality conjecture states that 
marginal densities (jl.6p in the scaling limit, when Aj = Ao + x% /n K (i = 1, . . . ,1) are universal 
(i.e. they do not depend on V). The scaling exponent k depends on the behavior of the 
equilibrium density p(X) in a small neighborhood of Ao- If p(Xq) / 0, then k = 1, if p(Ao) = 
and p(X) ~ I A — Ao| Q , then k = 1/(1 + a). 

For /3 = 2 universality of local eigenvalue statistics was proved in many cases. For example, 
in the bulk case (p(Ao) / 0) it was shown in [TO] (see also [H]) that for a general class of V 
(the second derivative of V is Lipshitz in some neighborhood of Ao) the scaled reproducing 
kernel converges uniformly to the sin-kernel. This result for the case of real analytic V was 
obtained also in [4]. Universality in the bulk for very general conditions on the potential V 
was proved also recently in [11]. Universality near the edge, i.e., the case when Ao is the edge 
point of the spectrum and p(X) ~ |A — Aq] 1 ^ 2 , as A ~ Ao, was studied in [4]. There are also 
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results on universality near the extreme point, where p(X) ~ (A — Ao) 2 , as A ~ Ao (see [3] for 
real analytic V and [15] for general V). 

The crucial difference between the case (3 = 2 and other j3 is that for j3 = 2 all correlation 
functions (|1.6p can be expressed in terms of the reproducing kernel of the system of normalized 
polynomials Pj = 7^ a;- 7 + . . ., (j = 0, . . . , n — 1) orthogonal on the real line with varying 
weight 

(A) := e~ nV{X) (1.10) 

/ pS n) (A)pl n) (A)u;W(A) dA = 8 j>k for j, fc > 0. (1.11) 

The orthogonal polynomial machinery, in particular, Christoffel-Darboux formula and Christof- 
fel function simplify considerably the studies of marginal densities (jl.6p . Moreover, asymp- 
totics of orthogonal polynomials Pn^ are known (see [4] for real analytic V and the 

recent paper [9] for non analytic V) and they can be used to prove bulk and edge universal- 
ity. 

For P = 1,4 the situation is more complicated. It was shown in [19] that the problem 
can be reduced to universality of some matrix kernels (see (jl.2ip . (|1.22|) below), which also 
can be expressed in terms of orthogonal polynomials (jl.lip . but to control their behavior one 
need to control the invertibility of some matrix (see Section [3] for more details). According 
to Widom [20], if the potential V is a rational function, then we need to control the inverse 
of some matrix of fixed size depending of V (e.g., if V is polynomial of degree 2m, then we 
should control some (2m — 1) x (2m — 1) matrix). Till now this technical problem was solved 
only in a few cases. In the papers [5JE] the case V(X) = A 2m (l + o(l)) (in our notations) was 
studied. Similar method was used in [7] to prove bulk and edge universality (including the 
case of hard edge) for the Laguerre type ensembles with monomial V. In [18] universality in 
the bulk and near the edges were studied for V being an even quatric polynomial. In [161 E] 
bulk and edge universality were studied for (3 = 1 and real analytic even V with one interval 
support a. 

But there is also a possibility to prove universality of local eigenvalue statistics by using 
another technique. In [18J Sojanovich made an important observation (see Remark 5 of [18] 
or Section 3 of the present paper) which allows one to replace the problem to control the 
Widom matrix by the problem to control E^ {n. _1 A/" ra [/]} for /3 = 1,2, 4. Thus the problem to 
study the correction terms of the order n~ l for E^ {n _1 A/" n [/]} becomes especially important. 

In a remarkable paper |10] Johansson studied the expectation and the variance of n~ l N n [f\ 
up to the terms 0{n~ 2 ). This allows him, in particular, to prove CLT for fluctuations of M n [f\. 
Unfortunately, his method works only in the case of one interval support a of the equilibrium 
density p and polynomial V with some additional assumption. 

In the present paper we generalize the idea of [10] to the case of real analytical V with 
one interval support of p, without any other assumptions. Moreover, we give a more simple 
proof of this result and apply it to the proof of bulk universality for f3 = 1, 4. 

Let us formulate our main conditions. 
Condition CI. The support a of the equilibrium measure density p consists of a single 
interval: a = [a,b], — oo < a <b < oo. 

Remark 1 It is easy to see that changing the variables M' = 2(M ^—I)/(b — a), in the 

case (i) we can always take the support a = [—2,2], 
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Condition C2. The equilibrium density p can be represented in the form 

p(\) = — P(\)QX 1 / 2 (X + iO), inf P(A)>0, (1.12) 

2lT Ae[-2,2] 

where 

X(z) = z 2 -4, (1.13) 

and we choose a branch of X 1 ^ 2 (z) such that X 1 ^ 2 (z) ~ z, as z — >■ +00. Moreover, the 
function v defined by |1. Tty attains its maximum if and only if X belongs to a. 

Condition C3. V is real analytic on a, i.e., there exists an open domain D C C such that 
a C D and V is an analytic function in D. 

Remark 2 It is known (see, e.g., f^j) that under conditions CI and C3 for any f3 the equi- 
librium density p of the ensemble \1.3\) has the form H.12\) - 111.13]) with P > 0. The analytic 
in D function P in (1.12]) can be represented in the form 

Hence, condition C2 states that P has no zeros in [—2,2]. Note also, that in the paper flUtf 
it was assumed additionally that V is a polynomial and P has no zeros on the real line. 

The first result of the paper is the theorem which allows us to control the expectation 
and the variance of linear eigenvalue statistics. 

Theorem 1 Under conditions CI - C3 for any analytic in D function f we have 

EpWn[f}} = J f(\)p(X)d\ 

1/2 \ 1 I f(z)dz I g'(C)d( , _ 2 



+ n\J~ l ) (2^)2 f C2d ^7% f Cd Pfc) + ^ nAf) , (1.15) 
where the contour is defined as 

Cd = {z : dist{z,a} = d}, (1-16) 
d is chosen sufficiently small to have all zeros of P(C) outside of Lid, 

«M-/«®£, (1.17) 

and r Ut p(f) satisfies the bound 

KAf)\<Cd sup 1/(2)1, 

2:dist{,2,cT}<2d 

with Cd depending only on d. 
Moreover, 

Bp {\M n [f} - Ep{Af n [f}}\ 2 } < C d sup \f{z)\ 2 . (1.18) 

z:dist{z,cr}<2d 
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One of the important applications of Theorem [T] (see discussion above) is the asymptotic 
of log Q n fi. Since the paper [2] it is known that 

n' 2 log Q n fi = -£ v + C(log n/n), 

where 

By = - J log — ^— p{X)p{p)dXdp - f V{X)p{X)dX. (1.19) 

But for many problems it is important to control the next terms of asymptotic expansion of 
log Q n ,0 (for applications see discussion in [8j, where the complete asymptotic expansion with 
respect to n~ l was constructed for the case /3 = 2 under assumption that V is a polynomial 
close in a certain sense to Vb(A) = A 2 /2.) 

Theorem 2 Under conditions CI - C3 for any /3 

n- 2 logQ ni/3 = n- 2 log Q { ^ + \pS V + \p (1-20) 
1 (, P\ 1 I (V(z)-z 2 /2)dz j; g' t (()d{ , ^__ 2 , 



+ n I 1 " 2 ) W f C2d XV^) f Cd PMz-C) + } ' 

where log Q^jg corresponds to the Gaussian case V = X 2 /2, £ v is defined by fO|) . §/? 
— ^/3£v , and 

1-t 



P t (X) = tP(X) +l-t, g t (z) = tg(z) + —{z - Vz 2 - 4). 
Remark 3 By the Selberg formula (see e.g. ^fM ) for the Gaussian case we have 

As it was mentioned above, Theorem [T] together with some asymptotic results of [4] for 
orthogonal polynomials can be used to prove universality of the local eigenvalue statistics of 
the matrix models (II. ip . We restrict our attention to the case when V is a polynomial of 
even degree 2m such that conditions C1-C3 are satisfied. Moreover we consider only even 
n. It is known (see [19]) that the question of universality is closely related to the large n 
behavior of certain matrix kernels 

K ntl (\, »):=(. ^' 1 \ X ^\ X , -f 5 7 l(A / } ) for /3 = l,n even, (1.21) 
\(eS nt i){X,n) - e(A - p) S n ,i\jJ>,X) J 

K nA (X,p) := (. S ^[ X f^ "f Sn>4(A,A*)\ for ^ = ^ ^ 



K (eS nA )(X,p) S nA (p,X) 

Here e(A) = ^sgn(A), where sgn denotes the standard signum function, and (eS ni/ g) (A, p) = 
L e(x — x')S n> p(x', y) dX' . Some formulae for the functions S n> p that appear in the definition 
of K n p will be introduced in (|3.3p . (|3.4|) below. In order to state our theorem we need some 
more notation. Define 

sin irt 



TTt ' 

J* " Koo(t) dt - e(£ - n) K^n - £) 
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Furthermore we denote for a 2 x 2 matrix A and A > 



\/X 1 \ „ (VX 



A (X) := VA U U 



-1 



yfX) \ \/A 

Theorem 3 Let V be a polynomial of degree 2m with positive leading coefficient and such 
that conditions C1-C2 are satisfied. Then we have for (even) n — > oo, Ao S 1 with p(Xo) > 0, 
and for f3 G {1,4} that 

V K nf ( A ° + f. A o + f ) = + 0(n-^), 

T*% { A ° + f < A ° + f ) = K ™^) + °^ 1/2 )' 

where q n = np(Xo). The error bound is uniform for bounded £, r/ and for Xq contained in 
some compact subset of (—2,2) (recall that supp p = [—2,2] by Condition CI). 

It is an immediate consequence of Theorem [3] that the corresponding rescaled /-point 
correlation functions 

Pt,l A + —,..., A + — , p ii4 A + —,..., A + — 

\ w v/n / \ m In 

converge for n (even) — )• oo to some limit that depends on f3 but not on the choice of V. 

The paper is organized as follows. In Section [2] we prove Theorems [T] and O In Section [3] 
we prove Theorem [3] modulo some bounds, which we obtain in Section |U And in Section [5] 
for the reader's convenience we give a version of the proof of a priory bound (|l,9p . 



2 Proof of Theorems [T], [2] 

Proof of Theorem [TJ Take n-independent e, small enough to provide that a £ C D, where 
a £ C M. means the e- neighborhood of a. It is known (see e.g. [13])) that if we replace in 
(jl.3p . (|1.5p and fjl .6|) the integration over R by the integration cr e , then the new marginal 
densities will differ from the initial ones by the terms 0{e~ nc ) with some c depending on e, 
but independent of n. Since for our purposes it is more convenient to consider the integration 
with respect to a £ , we assume from this moment that this replacement is made, so everywhere 
below the integration without limits means the integration over a e . 

Following the idea of |10j . we will study a little bit modified form of the joint eigenvalue 
distribution, than in (|1.3p . Namely, consider any real on a and analytic in D function h(() 
and denote 

v h (() = v(0 + h((). 

n 

Let p n ,i3,h, ^/3,h{- ■ ■ }, PiBh ^ e ^ ne distribution density, the expectation, and the marginal 
densities defined by (jl.3p . (jl.5p and (|1.6p with V replaced by Vh- 

By (11.3P the first marginal density can be represented in the form 

/n 
e -neV h (X)/2Y[ lx _ Xil e e -neV h (X l} /2 "Q \ Xi _ Xj f dX 2 . . . dX n . (2.1) 

i=2 2<i<j<n 
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Using the representation and integrating by parts, we obtain 

Here 0{e~ nc ) is the contribution of the integrated term. In fact all equations below should 
contain 0(e~ nc ), but in order to simplify formula below we omit it. 

Since the function p^l h (X, m) is symmetric with respect to X,u, we have 

p£U( A > f^dXd/j, _ f pQ A (X, fJ>)d\dn f P^l h (X, fi)dXdu _ f pQ jh (X, fi)dXdu 



(z - A)(A - fx) J (z-A)(A- / u) J (z-^-X) J {z — X)(z — n) 
Hence, equation (|2.2|) can be written in the form 

KWpHhW dX = 2_ I PilhW dX + (n-l) f P { A h (X,»)dXdn 



z — X (in J (z — A) 2 n J (z — X)(z — u) 

Let us introduce notations: 



[ pW^ridXdy 2 f f P ^l h (x)dx y f&M 



dX 



(2.4) 



fcn,/3,fe(A,/i)dA^ 

(z-A)(2-p) 
where 

k n ^ h {X,n) = n(n - l)pQ A (X,n) - n 2 p { ^ h (X)pQ h (fj) + n<5(A - ^)pQ jh (X). (2.5) 
Moreover, we denote 

9n,p,h{z)= , V(z,X) = . (2.6) 



z—X z—X 

Then equation (|2.2j) takes the form 

g 2 n ^ h (z)-V'(z)g M (z) + j V(z,X)p { $, h (X)dX 

_ 1 fh'(X)p^l h {X) 1/2 \ fP^Ux) 



n 



f h (X)p\^ h {X) i/ 2 \ f p\^ h (X) I 
/ z-X dX - n U " ) J J^W ~ n^ h{z) - (2 - ?) 



Using that V(z, £) is an analytic function of C in D, we obtain by the Cauchy theorem that 
for any z outside of Cd 

V(z,X)p { "i h (X)dX = ^-.i V(z,Og nAh (Qd(, 
Thus, (|2.7|) takes the form 
9lM z ) ~ V'{z)g nAh {z) + V(z, ()g n MOd( 

n 



1 /*'(A)," ,(A) 1/2 \ /-p n ii(A) i 
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Passing to the limit n — > oo, we obtain for any fixed z the quadratic equation 

g\z)-V'(z)g(z) + Q(z) = 0, Q{z) = ^ / U(z,CMCR, (2.9) 



where g is defined by (|1.17p . Hence, 



g{z) = \v\z) - ^V'{zY-AQ{z). 

Using the inverse Stieltjes transform and comparing with (|1.12p . we get that 

2g(z)-V(z) = P(z)X 1 / 2 (z). (2.10) 

where X(z) is defined by (jl . 13f) . 
Denote 

Un,l3,h( Z ) = n (9n,/3,h{ Z ) ~ 9(z)) & 9n,/3,h{ Z ) = 9{ z ) + ~Un^,h( Z )- i 2 - 11 ) 

Then, subtracting (I2.9P from (|2.8p and multiplying the result by n, we get 

(2g(z) - V\z))u nAh {z) + ^~<f V(z, Ou nAh (()d( = F(z), (2.12) 

where 

F(z] = J ^^M^ + (2 _ ^ ^ w + i_ u , n ^ (2 13) 

--ul,/3,h( z ) ~ ~5n,pA Z )- 
n ,y, n 

Using (12TTUD . we obtain from (12321 

P{z)X l / 2 {z)u n ^ h {z) + Q n {z) = F(z), Q n (z) = ^~f V(z, Qu 

Then, choosing d such that the contour defined by (|1.16p does not contain zeros of P(C), 
we get for any z outside of 

^- £ (p(Qx i i\Q UnAh (o + Gn(o - no) p ic) d ( c z _ c) = °- ( 2 - 15 ) 

Since, by definition (|2.14|) . Q n (C) is an analytic function in D, and z and all zeros of P are 
outside of C^, the Cauchy theorem yields 

Sn(CR _ Q _ 



2mJ Cd P(C)(*-C) 

Moreover, since 



VAftO) = ~ z (/ dX Pilh{>) - f d\p{X) \ + nO(z~ 2 ) = n0(z- 2 ), z -»• oo 

we have 

^(^(^nOfr 1 ). (2.16) 



Then the Cauchy theorem yields 



27ri J c 

Finally, we obtain from (|2.15p 



i / x^jcWMOdC Y i/2 ( s , v 



Un,P,h( z ) 



(*-C) 



F(Qd( 

2iriX 1 / 2 {z) J c . P(Q(z-Q- 



(2.17) 



Now take d small enough to have all zeros of P outside of C-^d- Then for any z : distjz, a} = 2d 
equation (I2.17P implies 



Un,f3,h( Z ) 



F(z) 



+ 



XV 2 {z)P{z) 2riXV 2 {z) fc 3d P(C)(C 



According to the result of [2] for any /? we have a priory bound 

Cn log n Cn 1 / 2 log n , Cn 1 / 2 log n 



dist {z, a} 



dist 6 {z,a} 



(2.18) 



2.19) 



dist {z, a} 

where C is an absolute constant. 
Denote 

M n (d) = sup \un,p,h( z )\ 

z:dlSt{z,a}>2d 

By (|2.16p and the maximum principle, there exists a point z : dist{z, a} = 2d such that 

M n (d) = \u n £ jh (z)\. 

Then, using (j2.18|) . the definition of F (see (12.13P ). and (12.19p . we obtain the inequality 

M n (d) < -CiMl{d) + C 2 logn, 
n 

where C\ and C 2 depend only on d, sup |P _1 (z)|, sup \n~ 1 h(z)\, and from C 

dist{z,cr}<3d dist{z,o-}<d/2 
of (I2.19p . Solving the above quadratic inequality, we get 



M n (d) > (2Ci)" 1 (n + y/n 2 -4CiC 2 nlog 



n) 



M n (d) < (2Ci)" 1 (n - yjn 2 -ACiC 2 n\ogn) 

Since the first inequality contradicts to (|2.19p . we conclude that the second inequality holds. 
Hence, we get 

sup \un,p,h( z )\ — 2C2logn + C(sup |/t'(A)| + dist~ 2 {z, a}). 

z:dist{z,cr}>2d A eo- e 



Note that the bound gives us that for any real analytic 93(C) 



n 



^A)(p£l fc (A)-p(A))dA 



1 

2vri 



(2.20) 



< w n sup \<p(z)\ + sup \h'(X)\ 
\zec 2 d ago- £ 



where 

w n = 2C 2 log n 

Now we are going to use the following lemma, which is an analog of Lemma 3.11 of |lUj . 
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Lemma 1 If \2. 20\) holds for any real h, and some 92 which is analytic in Di C D (a £ C 
then there exists an n-independent constant C* such that 



The lemma was proved in [10J, but for convenience of readers we give its proof at the end of 
the proof of Theorem [TJ 

Applying the lemma to (pi (A) = — A) and <pz(X) = ^s(z — A) with dist{z, a} > d, 
and using (|2.20j) . we obtain for such z (cf (12. 19h ) 

\SnAh\ < C' d log 2 n, \u nAh {z)l\u' nAh {z)\ < C' d \ogn. (2.22) 

Then, using this bound in (|2.18j) instead of (|2.19p . by the same way as above we get (12.20j) 
with w n = Ci(sup Ag(7e |^'(A)| + C d ). Then, applying Lemma [JJ once more, we obtain that 

\6 n ,p,h\ < C'L \< Ah {z)\ < C'l (2.23) 

Using these final bounds in (|2.17p . we obtain that 

where 

\r n (z)\ < rT x C d . 

□ 

Proof of Lemma [TJ Take any real analytic ip such that sup 26 £ 2d |y(^)| < 1- Using the 
method of [JO], consider the function 



F n (t) = Ep th jexp 
It is easy to see that 



t 

2w, 



n ,. 

-Y^^Xi)- \ v(X)p(X)dX) 

i=l 



^ log F n (t) = (2w n y 2 B^ h+t ^ /2Wn I (lL&( X i) ~ ^,h+t V /2w n W(Xi)})^ J > 0. (2.25) 
Hence, by (l230l) . for t G [— 1, 1] 

log i? n (i) = logF n (t)-logF n (0) = I ^- log F n {r)dT <\t\^- log F n (t) 

Jo dr dt 

= \t\{2w n )~ l V^ h+tv/2Wn j£ ^p(Xi) - J <p(X)p(X)dX^J j 

= £/^ (A) {Pil h+ t,/2 Wn W ~ PW) dX < \t\ 

Thus, for t € [-1,1] 

F n {t) <S< 3, 

and for any t £ C, \t\ < 1 

\F n {t)\ < F n (\t\) < 3. (2.26) 
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Then, we have by the Cauchy theorem, for \t\ < I 



\K(t)\ 



F n (t')dt' 



1 



2vr L,, =1 (if -if 



<6, 



and therefore for \t\ < A 

11 — 15 



|F n (t)|= f(o)- Tf; 

•/ 



> 



Hence, log F n (i) is an analytic function for \t\ < A and so, using the above bounds, we have 



|t|=l/12 



fog F TO (t) 
t 3 



< C. 



Finally, using (|2.25j) . we get 

□ 

Proof of Theorem [2] Consider the functions Vf of the form 

V t (X)=tV(\) + (l-t)Va(X), (2.27) 

where Vb(A) = A 2 /2. Let Q n ,p(t) be defined by (|l,4p with V replaced by Vt. Then evidently 
Qn,p(l) = Qn,/3 and Qn,/s(0) corresponds to the Gaussian case Vq(A) = A 2 /2. Hence, 



ti ft tl j q at 



f /"W dX(V(X)-V (X))pQ(X;t), 



(2.28) 



where p^AX^t) is the first marginal density corresponding to V%. Using (|1.7p one can check 
that if we consider the distribution (|1.3p with V replaced by Vt, then the limiting DOS pt 
has the form 

ft (A) = tp(X) + (1 - t)po(X) = iV4-A 2 [*P(A) + (1 - t)Fo(A)] , (2.29) 

with X defined by ()1.13|) and Pq(A) = 1. Hence, using (|l,15p for the last integral in (|2.28|) . 

we get 



1 



72" 



log Q n fi 



n 



\ log QnA°) ~ f £ + f fv 



1 



72. 



P 



1 



(V(z) - z 2 /2)dz 



2 J (2^) 2 Jr XVH 



+ 0(n 



-2\ 



where Ey is defined by (|1.19p . F t and gt are defined in (|1.20p and £q = — § (see, e.g., p2]). □ 
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3 Bulk universality for orthogonal and symplectic ensembles 

In a remarkable paper [19] Tracy and Widom showed how to express the functions S n p that 
appear in the definitions (jl.21j) . (jl.22j) in terms of orthogonal polynomials defined by (|l,10p - 
(jl.lip . Set ipj := PjVvffi, j > 0. Then the system {ipj}j>o defines an orthogonal basis 
in L 2 (M) with respect to the standard inner product (/, g) := f R f(X)g(X) dX. Moreover, they 
satisfy the recursion relations 

M U) (A) = (A) + b^t\x) + afrW&M, (3.1) 

which define a semi- infinite Jacobi matrix J^ n \ It is known (see, e.g. [13] ) that 

l4 n) |<C, \b^\<C, \n-k\<en. (3.2) 

(n) 

In order to state the formulae for S n> p we need to introduce more notation. Let D^o and 

(n) 

M<£ be semi-infinite matrices that correspond to the differentiation operator and to some 
integration operator respectively. 

AfW : = (e$ B) ,v£°) with (e/)(A) := / e(A d». 

Both matrices Doo and Mso are skew-symmetric. Using in addition that for j < k 
(D&l) jk = Jj( P f(X)y-^V\X) P f(X)) p P(X)w^(X)dX 

= -~l y\x) v f{x) v t\M n \x)dx = -^v\j^) ]k 

by orthogonality and the spectral theorem, we see 

(Dt ] )j,k = for \j -k\>2m and 
\(D^) jik \ < nC, \j-n\,\k-n\<en. 
In particular, we may write 



j,k>0 



)jkV k 

\k-j\<2m 

as a finite sum. Since )' = we have for any j,l > that 

% = ( e (^ n) y,^)= E (4 n) k(M< R) )«. 

|fe— i|<2m 

(n) fn) 

This relation together with the skew-symmetry of Moo and Doo proves 

£)WjVf( n ) = 1 = M (n) D (n) 

-^oo -'"oo JV - I oo -^oo ■ 

Next we denote by the principal n x n submatrices of and Z?oo\ i.e. 



M<»> := ((MW) ifc ) < i) ,< n _ 1 , DW := ((4 n) )^)o<, 



fc<n-l- 
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The formula of Tracy- Widom for S n 3 now reads 



n-l 

Sn,i(\n) = -E ^" ) (A)(M("))- fe 1 (6^ n) )(/i) (3.3) 

i,fc=0 
n-l 

Sn/ 2 , 4 (A,/x) = - £ (^(AX^jM^OO (3.4) 

j,k=0 



As a by product of the calculation in [19] one also obtains relations between the partition 

Qn/2,4^r 



functions Q n ,p and the determinants of and 



detMH = (^§) , det# 



(n/2)!2«/ 2 



where T n := Ilj=o 7/ an< ^ * s ^ ne lading coefficient of pj 1 ■ It is also known (see [12]) 
that Q n ,2 = r n / n ' an d we arrive at 

Since = 1 and [D^) jk = for |j - k\ > 2m - 1 we have = 1 + A n with 

A n being zero except for the bottom 2m — 1 rows. Define T n to be the (2m — 1) x (2m — 1) 

(n) (n) 

block in the bottom right corner of D n M ra , i.e. 

(T n ) jk := {D^M^) n _ 2m+j , n _ 2m+k , 1 < j, k < 2m - 1. 
lat de 

by Stojanovic in [18] : 



Then we have that det(T„) equals det(M n n ^ D^) and we arrive at a formula, first observed 



det(T n ) = ( ® n >) Q f^ ) • (3.5) 
V ! \Qn,2(n/2)\2n J V ) 

Since equals 1 up to the matrix A n of rank 2m — 1 (independent of n) it is 

conceivable that one may express (M^)" 1 and (D^)- 1 that appear m (035 by # 

and -Bn™^ respectively up to some correction terms that involves the inverse of T~ x . Using 
this idea Widom provided in [20J a useful formula for S n p that was later refined in [Tj- In 
order to present this formula introduce some more notation: 

fo( n ) — ( I ( n ) / ( n ) / ( n ) \T 

^1 :— vPn-2m+H V n -2m+2' • > ^n-1/ 1 ' 

$W — (./» ,/.(«) ,/.(") 

^2 ■- > V-n+1' • • • ' rn+2m-2J 

and 

M rs := (£*(»), ($i n) ) T ), L>, s := ((<^ n) )', ($i n) ) T ), 1 < r,a < 2 

define some (2m - 1) x (2m - 1) submatrices of and . Observe that = 1 

together with (-D00 )jjfc = for |j — k\ > 2m implies 

T n = 1-D 12 M 21 . 
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Finally we denote by K n (X, p) := Y^j=o V^MV^O-O * ne reproducing kernel. We then have 
® 

S n>1 (\,n) = K n (X,fj,) + <S> 1 (X) T D 12 e$ 2 (p) - ^(XfGe^ip), 
G := D 12 M 22 (1 - D 2X M X2 y x D 2X 
S n /2,A{\ M) = K n (X, p) + ^ 2 (X) T D 12 e^(p) - $ 2 (A) T Ge<f 2 (/i), 
G := -D 21 (l - M 12 D 21 )- 1 M 11 D 12 

Since S n> p(X, p) = — S nt g(p,X) one has for even n 

(eS n ^)(X, p) = - S n ,i(t,p)dt, {eS n/2A ){X, p) = - / S n/2A (t, fj) dt. 
Jx Jx 

Using this representation in the 21-entry of K n g together with detT n = det(l — D 2 \M\ 2 ) = 
det(l — Di 2 M 2 i) its straightforward to see that Theorem [3] follows from the following Lemma. 

Lemma 2 Given any compact set K C (a, b) there exists a C > such that for all n > 2m 
and all j, k S {n — 2m + 1, . . . , n + 2m — 2} one has 



(a) SU p|e^ n) (A)|<^; (b) \(M^) jk \<-; (c) | logdet(T n )| < C. 



C fm ^ , C 



xek' 3 ' ' \'" ' ' n 

Statements (a) and (b) will be derived from the asymptotics of the orthogonal polynomials 
in Appendix 01 We now prove statement (c) using Theorem 1. 

Consider the functions Vt defined in (|2.27|) Then, as it was mentioned above (see the 
proof of Theorem [2]) the limiting equilibrium density pt has the form (|2.29p . Hence, for any 
t £ [0, 1] Vt satisfies conditions C1-C3 and if we introduce the matrix T n {t) by the same way 
as above for the potential Vt, then T n (0) corresponds to the GOE and GSE. Consider the 
function 

L(t) = logdetT n (i). (3.6) 
To prove that |£(1)| < C it is enough to prove that 

|L(0)|<C, \L'(t)\<C, t€[0,l] (3.7) 

The first inequality here follows from the results of [19] . To prove the second inequality we 
use (I3.5P for V replaced by Vt. Then we get 



L'(t) = n 2 j AV(X)p^lf(X)dX + n 2 j AV(X)p { ^l t (X)dX - 2n 2 j Ay(A)^" 2 ) 4 (A)dA. 
It is easy to see that 

h > Phi? = lim Phlt = lini Phlt = PtW 



with pt defined by (I2.29p . Hence, using (I1.15p . we obtain that the first and the second terms 
of (|1.15p give zero contributions in L'(t), and therefore 

L'(t) = 2(r n/2At (AV) + r nXt (AV) - r at (AV)). 

But, according to Theorem [U all terms here are bounded uniformly in n. Thus, we have 
proved the second inequality in (|3,7p and so statement (c) of Lemma (2j 
□ 
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4 Appendix: uniform bounds for (M$)ij 

Set 

5 n = n - 2/3+K , < k < 1/3. (4.1) 

Then, according to [4J, we have 
«A) = + |A|<2-* n ; 

^(A) = C f^§ d + OCn- 1 )), |A|<2- 5n ; 

4™)(A) = n^P^Ai (±n 2 / 3 $ ± (A T 2)) (1 + 0(|A T 2|)) (4.2) 
+n- 1 / 6 B^ ) Ai' (±n 2 / 3 $±(A T 2)) (1 + 0(|A T 2|)) + O^- 1 ), |A T 2| < <5„; 



V&(A) = n^P^A* (±n 2 / 3 d> ± (A T 2)) (1 + C(|A =F 2|)) 

+n~ 1 / 6 pg ) Ai / ((±n 2 / 3 $±(A T 2)) (1 + 0(\X T 2|)) + 0(0, |A T 2| < 5 n ; 
|4")(A)| < e —(W- 2 ) 3/2 , |^ l(A) |< e -nc (|A|- 2 )3/ 2) | A | >2 + Jn . 
where 

1 /"^ 1 2 

p n (A) = - / P(A)\/4- A 2 (fA + - arccos(A/2), P n _i(A) = F„(A) - - arccos(A/2) (4.3) 

2 Jx n n 

with P defined in (]1.14p . Functions &± in (|4.2p are analytic in some neighborhood of and 
such that <3?±(A) = a±x + 0(x 2 ) with some positive a±. 
Denote 

sinnP n (A) sinnF n (2-5 n ) sinnF n (-2 + 5 n ) 

A \ X > ~ ~ 1 |A|<2-5 n H 1 |A-2|<5„ H 1 |A+2|<5 n 

n-V" (* (n 2 / 3 $ + (A - 2)) - * (n 2 / 3 $ + (-^))) 1, A _ 2 |< 5 „ 

_ n -i/2 (_n 2 / 3 d>_(A + 2)) - * (-n 2 / 3 d>_(0)) l|A+2]<5 n (4.4) 

1 pJ+ } Pj") 

5(A) = PA(A)XV4( A ) 1 W^- a » + $' + (A - 2) 1|A - 2 ^ a » + $'_(A + 2) 

:= I Ai(t)dt. (4.5) 

■/ — oo 



with 



Proposition 1 Under conditions of C1-C3 for any smooth function f we have uniformly 
in [-5 n - 2,5 n + 2] 

e(M n) )(A) = A(A)P(A)/(A) + er n {\) + O^ 1 ) + l, A±2 |< 5n O(n- 5 / 6 ), (4.6) 

where 

2+S n 

\r n (X)\dX < Cn- 1/2 - 3K/A . 

—2—8 n 

Similar representation is valid for ^(fip^i) if we replace in ^4-4\ ) F n by F n _i and B^' by 
Pgl • Moreover, it follows from that 

|6(M n) )(A)| < Cn- l '\ \e{f^l x )[X)\ < Cn- l '\ (4.7) 
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Proof. We use the following simple relation, valid for any continuous piecewise differentiable 
functions A, f, and any piecewise differentiable B, if A(X)B(\)f(X) — > 0, as A — > ±00: 

e(A'Bf)(X) = A(X)B(X)f(X) - e(A(Bf)')(X), 

where (Bf)' may contain 5-functions at the points of jumps of B. By the choice of A, B (cf 
(|4.4p and (14,21) ) e(A'Bf) corresponds to the principal part of e{fijj < if > ). The terms C(n _1 ) 
and l\\±2\<s n O(n~~ 5 / 6 ) in (|4.6p appear because of the integrals of C(n _1 ) in the second line 
of (|4.2p and the terms in the forth line of (|4.2p respectively. Hence we need only to prove the 
bound for r n = A(X)(Bf)' . Observe that 

r2+8 n r — 2+<5 n N 



JK(X)\dX < ^[J 9s +J_ 2 _ g J\(Bf)'(X)\dX 

1 r2—8 n 

+- / |(B/)'(A)|dA = n-^OiSn) + n" 1 ©^" 3 / 4 ) = ©(n" 1 / 2 " 3 ^ 4 ). 



2— 8„ J— 2— 8„ 
2—8 n 

□ 

Using recursion relations (|3.ip it is easy to get that for any \j\ < 2m 

iffi^A) = /q;(A)^(A) + fum^iW, 

where foj and /ij are polynomials of degree at most Note that since it is known that 
and for k — n = o(n) are bounded uniformly in n, foj and fij have coefficients, bounded 
uniformly in n. Hence for our purposes it is enough to estimate 

h := (Jl3^2i,e(J 0k ^)), h := ^3 := (/oi^ n) , e(/ rf>)). (4.8) 

It follows from Proposition Q] that 

A = hi + h2 + hs + er n ) + ©(n- 1 ), 

where 

_! f 2 -^ / 0J - (A)/ifc(A) sin nF n (A) cos nF n _i(A) 
711 = n ^(A)(4-*)V ^' 

J 12 = ^WfiW (n 2 / 3 d> + (A - 2)) - * (n^ + (-5 n ))) 

Ai(n 2 / 3 $ + (X-2)) , s , 
" $V(A-2) - ^W/ifc(A)dA, 

and I13 is the integral similar to I12 for the region |A + 2| < 5 n . It is easy to see that 

j P (+) p (+) /oj( 2 )/ifc( 2 ) n 1 , n u 

/l2 " Sl1 521 2n(d>' + (0)) 2 (1 + ' 



I 
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Sl1 ^ 21 2n(<I>' + (0)) 2 + 



Moreover, using the bound for r n from (j4.6|) and (|4.7p . we get 

(/U^i.cr,,) = -(c(/y^i).r») < / K/y^i)l |r„| dA = ©(n" 1 - 3 */ 4 ). 
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Hence we are left to find the bound for In. 

T (0 f 2 - 5n f 0j (X)h k (X)smn(F n (X) - F n ^(X)) 

hl = (2n) Us, ^n(A)(4-A')V2 dX 

, f 2 - S " foj(X)f lk (X)smn(F n (X)+F n ^(X)) ^ _ 

+( n) L +Sn Kjx^>?W dA - J n+% 



By the definition of F n and F n _i (I4.3p . we obtain 

f-2 



= i + (i) r 2 f 0j (x)f lk (x) 

11 2n 7_ 2 P(A)(4-A 2 ) 1 



/2 dX. 



Moreover, integrating by parts one can get easily that I'd = 0(n 2 S n 3 ^ 2 ) = 0(n 1 3k / 2 ). 
The other two integrals from (j4.8j) can be estimated similarly. 
□ 

5 Appendix: proof of the bounds (12.191) 

Let us introduce a function H which we call Hamiltonian to stress the analogy with statistical 
mechanics. 

n 

H(A) = -Y,(V(X l ) + n- 1 h(X l )) + 2 £ log | A, - A,|, A = (Ai, . . . A n ). 

i=l l<i<j<n 

It is evident that for any continuous /(Ai, . . . ,X k ) 



/"m x W1 , x / /(^ ■ ■ ■ | X k )e^ H ^dA 
J J(Xi,...,* k )p kJ3:h (M,---,X k )dA 1 ...ciA k - J e n/3H(A) dA _: 



Moreover we introduce the "approximating" Hamiltonian, depending on a functional param- 
eter m : supp m C [—2, 2] 



H a (A;m) = ^\ n (Aj;m) + (n - l)£[m,m]. 



1 n — 1 

v n (A;m) = -V(X)--h(X) + 2 £(A;m), (5.1) 

n n 



Here 



£(A;m) = J log |A — fi\m(fj,)dfj,, 

C[m,m] = J log | A — fj,\^ 1 m(X)m(n)dXdfi. 

By the Jensen inequality for any two real functions 7i\(A), ^(A) we have 

Je^^dA enmm -n^ Je^^dA ^-H^ 

J e n/3H2(A)/2 dA ~ ' f e nPHl(A)/2dA ~ 

where we denote (. . . ) m& =/(... )e n/3n ^ A ^ 2 dA/ f e npH *^ 2 dA (5 = 1,2). Then we get 
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Taking here H\ = H, H2 = H a , we obtain 

_ J(H a - H)e-P nH / 2 dA J{H a - ff) e -/^(A;m)/2 dA ^ 
R[ml : " (n-l)/e-^/2dA " (n-l)/e-^(A; m )/2 dA " : ^M' ^ 2 > 

Since H and i7 a are symmetric, we can rewrite the l.h.s. of (|5.2p as 

m) - pS^WpS,^)) dAd P + ^bf/U - m ^i?/3,h - H 

(5.3) 



i?[m] = / logjy 



where p^ n ] h and ] h are defined by (|l,6p if we replace V by V^. To obtain the expression 
for the r.h.s. of (|5.2j) we need to replace p^{X) and p^(X, p) in (]5,3p by pj"g a ^(A;m) and 



Pi /3 /A A ' m )Pi ,3 /KM! m )> ~~ correlation functions of the approximating Hamiltonian ([5]), where 
P { iSl^ m ) = e^ A;m >/ 2 ( J dXe^^ 2 ^ . (5.4) 



This yields: 

Now let us choose the function m. Set 

n 

n — 1 \ r v y ' /3n 



JittH = ^[PiXl - m '^i?/3?fc - m L ( 5 - 5 ) 



m «( A ) = I^T (pW + ^n(X) ) 1|A|<2, (5-6) 



where 



^n(A) 



^ ^ mKM + m.)m + _ ,, 1)(A) + a „,,», (A)i (5 . 7) 

^ ./-2 (/i-A)V4-/i 2 7rV4-A 2 



the function (log p) n (X) coincides with log p(A) on the interval <r n = [—2 + n~ 1//2 , 2 — n^ 1 / 2 ], 
and (logp) n (A) is a linear function for A 6 a \ cr n , chosen so that (logp) n (A) has continuous 
derivative on a. The constant a n here is chosen to provide the condition 



J m n (X)d\ = l^a n = -0 - J vjp (X)dX 



B _ f 2 ((logpyM + ph'irifflndy = 1/4 

7-2 v 7 !^! 2 " 1 ' 



Since p has the form (11 . 12H . iffi is a sum of a bounded function which comes from P and of 
a negative function which comes from the integral of (log\/4 — A 2 )'. Hence, 

|i/W(A)|dA <C- f vjp (A)dA = ©(n 1 / 4 ). (5.8) 



It is easy to see that ^ n (A) is chosen to satisfy the equation 

2 v n {p)dp 



-2 A ~P 



(logp)UA)+/3/ i '(A)/2. 



Therefore 

£(A, i/ n ) = (log p) n (A) + Ph(\)/2 + r n (A), (5.9) 
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where for |A| < 2 r n (A) = C n and C n is a constant independent of A, but depending on n. 
One can find C n as 

C n = £(0,O - (logp) n (0) - (3h(0)/2 = J log|AK(A)dA - (logp) n (0) - /3/i(0)/2 = Oin 1 ^) 

(5.10) 

(here we used that i>n for |A| < 1 is bounded uniformly in n ). Hence, 

(5.11) 



f3nv n (\,m n )/2 = Pn[£(\,p)-V(\) /2 + (logp) n (A) + C„, |A| < 2 



Let us estimate ^£(A,^ 1} ) for A > 2. From ([521) we get 



dA 



^(AiJdAi 



A-Ai 

2 



v / 4^Af((logp)U^) + /3/ i , (/i)/2) 



l r r- 

— / dX\ / dit- 

W-a 7-2 (A-Ai)(m- Ai)^4^7? 



l v / A^y(iogp);(») + /3//( M )/2) i 



< su P {|(logp)» + pti(ji)/2\} f dAtfl + ) (4 - n 2 )- 1/2 d^ < n^d. 



A — /i 



Here we used the identities (valid for A ^ [—2, 2] 



V4^Af 



V / A 2 ^~4 1 f 2 (4-u 2 )- 1 ^ i 



A — /i 



^7-2 (A - Ai)(/i - Ai) A-/i 7Tj_ 2 

and the bound sup |(logp)^(/i)| < Cn 1 / 4 . Therefore 

£(A, i/ n ) - C{2, u n ) < dn l / A \\ - 2| + Can 1 / 4 ^ - 2] 1 / 2 , 

where the second term in the l.h.s. comes from a n v^ in (|5.6p . Moreover, since under 
conditions C2, C3 there exists C* such that 

£(A, p) - F(A) = C*, |A| < 2, £(A, p) - F(A) - C* < -C \X 2 - 4| 3 / 2 , | A| > 2, 

we obtain 

n(3v n (\,m n )/2 - n(3C*/2 - C n = (log p)„l| A |<2 + r n (A)l| A |>2, 
r n (X) < Cn 1 / 4 (|A| - 2) 1 ' 2 - nC \X 2 - 4| 3 / 2 . 



Then we have 



nPv n (\,m n )/2-nf5C* /2-C ndx < 2 / e Cn ^ C ° nx3 d,X < (J n ~2/S 

Jo 



(5.12) 



The last bound can be obtained by splitting the interval [0, oo) in two parts: [0,re _2//3 ) and 
[n~ 2 / 3 ,oo). Then in the first interval we used the fact that sup{Cn 1//4 x 1//2 — nx 3 ^ 2 } < c 
and in the second interval we used that this function is negative, its derivative is a negative 
decreasing function, bounded from above by (— Cn 2 / 3 ). 
For |A| < 2, since 

| e (iogp)„ _ p \ < C n -i/4i CTW we g et 



(logp)n dA = / p ( A ) dA + o(n" 3 / 4 ) = 1 + 0(n" 3 / 4 ). 
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Hence, using the above inequality and (|5,12p . we obtain 

Thus 
and 

£\Plj3%(^ m n) ~ m n,Pi"^i(A; m n ) - m n ] < C (n~ 2 C[v n , u n ] + n~ 4/3 + n~ 3/2 log; 
Moreover, using (|5.9|) and (|5.8jl . we write 

L[u n ,u n ] = -J {[\og P ) n (X) + ^h{X)/2 + C^jv n {X)dX<Cn 1 / 2 . (5.13) 

Finally we get 

^iP^ a J l (^' m n) ~ m n ,p^ h {X\m n ) - m n ) < Cn~ 4/3 . 
The inequality combined with (15. 2p gives us 



log |X^(^U A >^ (5-14) 
Let us prove that 

log P^( P ?^ (A '^ -Pi?i,&( A )Pil/»(^) dAd A* > -Clogn/n. (5.15) 

Introduce the function 

l n (X) = log | A|- 1 l| A | > „-r + (logn 7 + n 7 (n- 7 - |A|)) l| A | <n -7. 

It is easy to check that for any / the Fourier transform l(k) > and hence, for any 
positive operator K : L2(IR) — > and such that J K(X, n)dXdfi = 0, we have 

J l n (X - n)K(X, P )dXdn > (5.16) 

From (jl.6p it is easy to obtain that for any |x| < n~ 3 

IpS3,a( A + *) " < Cn-Hsup \V'\ + sup |/»'|), 

therefore 

1> I p^l h (X)dX>(l-Cn-\sup\V'\ + n-Uup\h'\))maxp { ^ h (X)n- 3 . 



Hence 
Similarly 

» 



maxpj n j h (X) < n 3 (l + Cn _1 (sup \V'\ + n" 1 sup |/»'|)). 



maxp [ ^ h (X,n) < n 6 (1 + Cn' 1 (sup \V'\ + n" 1 sup |/i'|)) 
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Then, the above bounds and the inequality 

/ | log tt— — c - l n (X - p)\d\ < Cn~ 7 log n, 
J |A — yu| 

imply 



log p^Tv^ft^^) -^ij,h( A Hj,h(^) J dXd P 

= ? 1 rr / Z„(A -//)*:„(*, AOdAd/x-^-r / / n (0)^ n > (A)dA - ©(n" 1 logn), 
n(n — 1)7 n — 1 J 

where the kernel A: n is defined by (12. 5p . Since A: n is positively definite, and J k(X, p)d\dp = 0, 
we can use ()5.16p . and taking into account that / n (0) = O(logn), obtain (|5.15p . 
Then 

\ 2 



And since it follows from (|5.6|) and (|5.13p that 

£[m n - p,m n - p] < Cn^^ 3 , 

we have 



ttfxlh ~ P'Pllh -P\< Cn~ l logn. (5.17) 

For any 3z ^ y, taking the Fourier transforms p^l h and p of the functions pf^ h and p, we 
get 

/ ^fl~ P(A) ^| < 2 |^(fc) - p(k)\e- Mvl dk 

< 2( / Ifcle-WW*) ^ ( / ^W-W 2 ^ = 2M-£V^ _ ^ _ p] . 

Then (|5.17p yields the second and the third bounds of (|2.19p . The first bound follows from 
Lemma [TJ 
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